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In this paper we consider the cosmological implications of dark energy models with a coupled 
system of a dynamical scalar field (the quintessence) and the neutrinos. By detailed numerical 
calculations we study the various possibilities on the evolution and the fates of the universe in this 
class of models. Our results show that due to the interaction with quintessence, neutrinos could be 
dominant over the quintessence in the future universe, however would eventually decay away. 



I. INTRODUCTION 

Recent data from type la supernova (SN la) JJ and 
cosmic microwave background (CMB) radiation 0] have 
provided strong evidences for a spatially flat and ac- 
celerating universe at the present time. In the context 
of Friedmann- Robertson- Walker cosmology, this acceler- 
ation is attributed to the domination of a new matter 
component with negative pressure, dubbed dark energy 
0. The simplest candidate for dark energy seems to be 
a remnant small cosmological constant or vacuum energy 
with p ^ {2 X 10^-^ eV)'^. This energy scale ~ 10"'^ eV is 
smaller than the known energy scales in particle physics 
except that of the neutrino masses, which is comparable 
to the scale of dark energy. 

The dark energy could also be due to a dynamical 
component, such as a canonical scalar field named 
quintessence P, |H IS- Cosmological observations indi- 
cate that the potential of the quintessence field should 
be very flat around the present epoch. Consequently the 
effective mass should be extremely small, mg ^ lO"'^'^ 
eV, which surprisingly is also connected to the neutrino 
masses via a see-saw formula mg ~ m'^/mpi, with rupi 
being the Planck mass. 

Are there any connections between the neutrinos and 
dark energy? Given the arguments above it is quite in- 
teresting to make such a speculation on this connection. 
If yes, in terms of the language of particle physics, it 
requires the existence of new dynamics and new interac- 
tions between the neutrino and the dark energy sector. 
Recently there are quite a few studies in the literature on 
the possible realization of the models connectin g ne utri- 
nos and dark energy^ 0, H H [H [H IH El lUlll El ■ 



Qualitatively these models have made at least two inter- 
esting predictions: 1) neutrino masses are not constant, 
but vary during the evolution of the universe; 2) CPT 
is violated in the neutrino sector due to the CPT violat- 
ing Ether during the evolution of the quintessence scalar 
fieldE3l- Quantitatively these predictions will depend on 
the dynamics governing the coupled system of the neu- 
trinos and dark energy. 

One of the possible couplings between the neutrinos 
and the scalar field is the derivative interactionllil: 



A 



(1) 



During the evolution of a homogeneous quintessence 
scalar field, 4> does not vanish and gives rise to CPT 
violation in the neutrino sector. However, since 4> is very 
small at the present epoch this type of cosmological CPT 
violation is predicted to be much smaller than the cur- 
rent experimental limits. But in the early universe with 
high temperature it has been shown in Ref . E3| that this 
CPT violation is large enough for the generation of the 
baryon number asymmetry via leptogenesis. This new 
mechanism ^^^^^for baryogenesis/leptogenesis provides 
a unified picture for dark energy and baryon matter of 
the universe. 

Another type of the interaction between the neutrinos 
and dark energy is that the quintessence field couples to 
the neutrino mass term. In the minimal extension of the 
standard model of particle physics, the neutrino masses 
can be described by a dimension-5 operator 



= jIlIlHH 



h.c. 



(2) 



where / is the scale of new physics beyond the standard 
model which generates the B — L violation, IltH are 



^ In analog to the idea of top quark condensate as a mechanism 
of generating the electroweak scale, one might think that the 
neutrino condensate gives rise to the scale of dark energy. In this 
sense the dark energy scalar field, such as quintessence behaves 
like a bound state of the neutrinos. Effectively the system of 



dark energy sector consists of neutrinos and (single or multi) 
scalars which interact with each other given (approximately) by 
an effective Lagrangian similar to Eqs. JHJ and (1). 
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the left-handed lepton and Higgs doublets respectively. 
When the Higgs field gets a vacuum expectation value 
< H >^ V, the left-handed neutrino receives a Majorana 
mass nil, ~ In Ref. |^ we considered an interaction 
between the neutrinos and the quintessence 



iLhHH + h.c, 



(3) 



where /3 is the coefficient characterizing the strength of 
the interaction between quintessence and the neutrinos. 
In this scenario the neutrino masses vary during the evo- 
lution of the universe and the neutrino mass limits im- 
posed by the baryogenesis are modified. 

The operator Q is not renormalizable, which in prin- 
ciple can be generated by integrating out the heavy par- 
ticles. For example, in the model of the minimal see-saw 
mechanism [13 for the neutrino masses, we have 

-Cncutrino = h.jluNRjH + ^M.jN'^^Nrj + h.C. , (4) 

where Mij is the mass matrix of the right-handed neu- 
trinos and the Dirac mass of neutrinos is given hy mjj = 
hij < H >. Integrating out the heavy right-handed 
neutrinos will generate the operator in ((2Jl. However, 
as pointed out in Ref. to make the light neutrino 
masses vary there are various possibilities, such as by cou- 
pling the quintessence field to the Dirac masses and/or 
the Majorana masses of the right-handed neutrinos. In 
Ref. llC| we have proposed a model of mass varying right- 
handed neutrinos. In this model the right-handed neu- 
trino masses Mi are assumed to be a function of the 

quintessence field Mi{(f>) = MiC^ "p' . Integrating out the 
right-handed neutrinos will generate a dimension-5 oper- 
ator like in (3), with the light neutrino masses varying in 
the following way 



Zfi^e ^'tijl^l^HH + h.c. 



(5) 



In this paper we will give a systematic study on the 
effects for a coupled system of light neutrinos and the 
quintessence field, especially on the late-time evolution 
and the fate of the universe. We will show numerically 
that 1) unlike the case in the absence of the interaction 
the relic neutrino energy density may not decrease, in- 
stead it can track the energy density of the quintessence 
field. In the future, all of the pressureless matter may 
be diluted away and only the components of the neu- 
trinos and quintessence remain. 2) Neutrino masses will 
increase as the universe expands, however when kinemat- 
ically allowed the neutrinos will decay for instance into 
electrons and pion mesons. Since the electrons and the 
mesons are not coupled to the dark sector, the fate of 
the universe is governed by the "pure" dynamics of the 
quintessence field. In this paper we will quantify the 
moment when the neutrinos decay by the numerical cal- 
culation. 



The paper is organized as follows: in Section II we 
will study the general features of the dark energy models 
with mass varying neutrinos and discuss analytically the 
behavior of the evolution of the system. In Section HI 
we present the numerical results for some specific models. 
Section IV gives discussions and conclusions. 



II. DARK ENERGY MODELS WITH MASS 
VARYING NEUTRINOS 

In this section we will study in detail the dark energy 
models where the neutrinos and the quintessence interact 
with each other. Specifically we start with a class of 
models given by the following Lagrangian 



(6) 



where Ci, — Di^v and £^ — ^df_,(j)d^(t) ~ ^{4') with 
V{(t)) being the quintessence potential; M{(t)) being the 
(/>-dependent neutrino mass such as in Eqs. ^ and |(SJl. 
In this paper we consider mainly the models of neutrino 
coupling to the quintessence via the mass term. 

As shown in Eq.(6) different dark energy models are 
specified by the form of and M{(p). In this paper 

we will consider two classes of models, one has the ex- 
ponential form of V{(l)) and M{(j)) and the other has the 
power law form. In this section we will give a general 
analysis of the behavior of the system and its condition 
for the attractor solution. 

The equation of motion for the scalar field is given 

by 



where 



dVi _ dM IM 
~d4 ~ ~d4^\E 



(7) 



(8) 



is the source term by the interaction between the neutri- 
nos and the scalar field, with n and E being the num- 
ber density and energy of the neutrinos respectively, () 
indicates the thermal average. For relativistic neutri- 
nos, the term ^ is greatly suppressed and the neu- 
trinos decouple from quintessence. For non-relativistic 
neutrinos, the effective potential of the system is given 
byT4//(0) = F(0)+nA/(<^). 

For this coupled system, the energy for each compo- 
nent does not conserve. It is easy to get that the energy 
density of neutrinos evolutes as 



(9) 



where M' — dM /dcf). Correspondingly due to the conser- 
vation of energy-momentum tensor of the whole system, 
the fiuid equation of the scalar field (j) is 



+ 3Hp^{l + w^) = -nM'4) 



(10) 
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where wj, — — = ^j, Y},'^} is the equation of state of 
the scalar field. We can define the effective equation of 
states w?/-f and w':^^ as 



+ 3Hp,{l + wlff) = , 



(11) 
(12) 



where 



and 



nM'^ _ l ^IogM 
" iHp, " 3 90 ^^^^ 



e// , l91ogA// 

^ 3 ad) Pd, 



(14) 



with = = -|. 

^ a log a /i 

The effective equation of states wf,^-^ and w'^^'^ describe 
the actual rate of the energy density decrease of the two 
components respectively as the universe expands. When 
the two equation of states become equal we expect the 
two components to remain a constant ratio in energy den- 
sities. In the future when the energy densities of all other 
components of the universe dilute away and requiring 



w 



eff 



w'l/f ^ we have 



(15) 



l91ogM 

where we have set SttG = 1 and used the relation = 

^ + r2y. We can use the requirement of a constant ratio 
of the energy densities to derive the attractor solution of 
the system. 

We first consider the class of models with exponential 
forms of V{(t)) and M {(()), which have an attractor solu- 
tion and will be shown below. For M{4)) — Me~^'^ we 
have ^ = — A. Requiring 



we get 



dloga 



(2(/)' - \)(t)" = 



(16) 



(17) 



This leads to the solutions (j)' = A/2 or 0" 0. The 
first solution gives Q.^ < 0. For 0" — 0, we have (j) — 
(j>o + A' log a. This solution leads to 



= -{\'^ - XX') 



(18) 



For the exponential potential V = Voe^'^, we have the 
equation of motion|l9j 

Y + (P™+3P.+2V^)- = -l3V+Xp^ . 

(19) 



From Eq. (|19|) we can get 



f +/3 + A 



Eqs. l(TS|l and (|^ give the unique solution: 

-3 



A' = 



X + (3 



(20) 



(21) 



This is the attractor solution given in Ref. [13 in the 
discussions of a coupled system of quintessence and dark 
matter. 

The exponential potential has very good attractor be- 
havior and generically the coupled system is very close to 
the attractor region today and leads to a constant ratio 
of the two coupled components' energy densities. This 
scenario seems to be less interesting to us, since the evo- 
lution of the universe is still dominated by the dynamics 
of the scalar field with a negligible neutrino component. 

In the following we will consider a class of power law 
models where M = M(j)-I^ and V{4>) = Vb0" , /3 and 
a are the two model parameters to be specified for the 
numerical calculations in the next section. In this class 
of models, = and we have 



3 30 



(22) 



Demanding that remains a constant during the evo- 
lution of the universe, one gets 



(23) 



In the next section we will present the detailed numerical 
results on the evolutions of the universe. Specifically we 
will focus on the power law models and on the behavior 
of the solution to the equation above, when tends to 
zero and (j)' /(j) tends to a constant. Our results will show 
explicitly the importance of the neutrino component and 
the interaction with quintessence in determining the evo- 
lution and the fate of the universe. 



III. NUMERICAL RESULTS 

In this section, with the power law forms of interaction 
and potential, we will present our numerical results for 
the evolution of the coupled system given in Eq. © . We 
will show that this class of models can satisfy the obser- 
vations within the model parameter space, furthermore 
predict different fates of the universe. At high tempera- 
tures, the neutrinos are relativistic and decoupled from 
the scalar field, as shown in Eq. ©. At low tempera- 
tures, the neutrinos begin to affect the evolutions of the 
scalar field. Specifically we notice that the system will 
undergo a period of oscillation before entering the at- 
tractor region. In the period of oscillation, the equation 
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FIG. 1: The f2i = ^ and efTective equation of state W as 
functions of the scale factor log a for Y = Vo(/>^ and M — 
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FIG. 2; The Qi and the effective equation of state W as 
functions of the scale factor log a for V = Vocj)'^ and M = 



of state of neutrino can be below —1, so its energy den- 
sity can grow up and become comparable with the scalar 
field or dark matter. Therefore, the fate of the universe 
as we show below can be quite different from the dark 
energy model in the absence of the interaction with neu- 
trinos. When neutrinos do not interact with dark energy, 
the density ratio between dark matter and neutrinos re- 
mains constant and typically dark energy will dominant 
the late time universe. However as shown below, the uni- 
verse can be dominant with dark energy and neutrinos 
after introducing the interaction, with constant ratio of 
the energy densities of the two components. On the other 
hand, the universe can also be dominated by neutrinos 
and dark matter with negligible dark energy in another 
scenario. In the numerical study we have taken the pa- 
rameter values today as h ^ 0.7, il^ « 0.7, r2,„ « 0.3 
and < 0.02. 

In Fig. m by specifying the potential V = Vo0^ and 
the neutrino mass M = M0~^^, we plot the evolution of 
Cli = y- and the effective equation of states as function 
of the scale factor logo. The effective equation of states 
in the figure are defined in Eqs. H13|l and H14|) . w total — 
Ptotal / Ptotai here is the equation of state of the whole 
system. After a short period of oscillation, the energy 
density fraction of the scalar field decreases while the 
component of neutrino increases to around 30%. The 
universe will be dominated by the two components in 
the future and the ratio of the two components keeps 
constant for ever. This picture corresponds to a solution 
of the forementioned equation H23|) when (p approaches 
zero and 4>' /4> remains nonvanishing.^ 

In the case of interacting neutrino and dark energy 



^ For simplicity we have taken neutrinos to be matter like in the 
early epoch. This does not change our late time picture when 
considering the neutrino behavior at the early times. 
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FIG. 3: The Q, and the effective equation of state W as 
functions of the scale factor log a for V — Vb0^ and M = 



there is typically a period of "build-in" phase before the 
system arrives the late time stable states. Fig. ^ stands 
for the case where the system enters quickly the stable 
phase. We find that to achieve a universe with a larger 
fraction of neutrinos, it will typically take a considerably 
longer time for the "build-in" phase. As an example we 
consider a model with V = VqcJ)'^ and the neutrino mass 
M = M(j)~'^. In Fig. [21we plot the evolution oi^li and the 
effective equation of states for this model. One can see 
after the period of oscillation, the universe approaches to 
a state consisting equally of the neutrino and the scalar 
field: — Q.tp — 50%. The equation of state of the 
system at the tracking region is w total = — 0.5. A more 
interesting case is plotted in Fig. 13 for V — Vbt/)"* and 
M = M(/)~^. The universe will finally be dominated 
by the neutrinos with U,^, « 2/3 while il^ « 1/3. The 
equation of state of the system is Wtotai ~ —1/3. 

The oscillating phase is determined by two conditions. 
One is the asymptotic behavior of V^^j:. In the case of 
a power law quintessence potential and neutrino masses. 
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FIG. 4: Wif, and the ratio of 1/(0) with neutrino energy density 
as functions of the scale factor log a for V = Votj}"^ and M — 



y^ff gets asymptotically to zero during the evolution. 
The other condition is the magnitude of V" compared 
with the Hubble expansion rate. The scalar field starts 
to oscillate at the time Vj^j ~ i?^. Due to the Hubble 
friction the oscillation damps with time. When the oscil- 
lating amplitude is small enough the kinetic term of cj) be- 
comes negligible, i.e. <^ ^^^d tends to — 1. 
For the power law form we will then have aV{(j)) ~ Ppu^ 

Taking the model in Fig. |31as an example, we delineate 
in Fig. ^the corresponding and V{4>)/ pi, as functions 
of the scale factor logo. The contribution of the kinetic 
term oscillates and decreases with time. After the freez- 
ing of W0 = —1, the ratio V{(j))/p„ approaches 1/2 and 
Vg'^j 0. Comparing with Fig. O we notice that wtotai 
and fl^ become nearly constant after approaches -1. 
However, we should point out that the damped oscillation 
lasts for a longer time and fi'^/il^ is not yet negligible for 
log a > 10 in the shown range of log a in Fig. Inland Fig. 
0] Meanwhile the two effective equations of state remain 
damped oscillating around Wtotai- 

As shown above in the class of models with the scalar 
potential V{(j)) and the neutrino mass M{(j)) being the 
power law form the neutrino masses keep increasing and, 
unlike the dark matter, its energy density does not di- 
lute away. In some cases neutrinos will even dominate 
over the quintessence. These results, however depend on 
the couplings of the neutrinos to quintessence. As an 
example we take ¥{(/)) = Vb(/)" while M — M exp(— /?</)). 
The behavior of neutrinos and dark energy cannot remain 
symmetric and their ratio cannot take as constant dur- 
ing late time evolution. The dark energy will transfer its 
energy density to neutrinos before their decoupling. The 
universe becomes neutrino dominant with a constant ra- 
tio between ft^ and flm- In Fig. |5lwe realize this picture 
with = Vo(t>^ and M = M exp{-4(j)/Mpi). The uni- 
verse evolves approaching a state dominated by matter 
and neutrino, where the effective equation of state Wtotai 
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FIG. 5: The Q and the effective equation of state W as 
functions of the scale factor log a for V — Vocj)'^ and A4 = 



V=V/,g=f' 

— V=V„t*,g=f 
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FIG. 6: The neutrino mass as functions of the scale factor 
log a, where g — M{(j))/M. 

is zero. Due to the different behavior of V/jy from the 
previous power law cases, the neutrino mass in Fig. [Slwill 
become constant in the future when the universe stops 
acceleration. 

Generally speaking the mass of neutrinos will vary as 
the evolution of the coupled system. It is easy to show 
that in general M(<^) = Ma~^'^'"' with a being the scale 
factor. For w'l/f < 0, the neutrinos are lighter in the 
past, but become heavier in the future. In Fig. El we 
plot the neutrino mass as a function of the scale factor 
logo for the models discussed in this section. It should 
be noticed that the neutrino mass will increase forever 
as the expansion of the universe for the power law cases. 
When the neutrinos are heavy enough, however they will 
decay to electron and pion ^ + . The life time 
of the decay is given by 

= 6. X 10-9 sec • (^1 - . (24) 

Therefore, the neutrinos will decay away immediately 
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once they cross the energy threshold. The energy den- 
sity of neutrinos will then be transferred to matter, which 
may eventually be diluted away with the expansion of the 
universe. In this case the scalar field will finally dominate 
the universe and determine the evolution of the universe. 
In Fig. El for an intuitive view, the moment when the 
decay occurs is set as — m^r + rrie- 

IV. DISCUSSIONS AND CONCLUSIONS 

In this paper we have studied the dynamics of the cou- 
pled system described by Eq. (|SJ). We paid particular 
attention to the effects of the back reaction of the ther- 
mal neutrino bath in the universe to the evolution of the 
scalar field — the quintessence — due to the i/i-dependent 
mass term'20j. 

The coupled system has many interesting phe- 
nomenologies. First of all, during the evolution of the 
universe the neutrino masses vary. This effect will make 
the matter power spectrum observed from the Large 
Scale Structure(LSS) survey different from that for a con- 
stant neutrino mass. Consequently the present cosmolog- 
ical limits on the neutrino mass rrii < 0.6eV will be re- 
laxed in the scenario with mass varying neutrinos 0| . In 
the current paper we mainly focus on the models where 
neutrinos can track dark energy. Although the mass of 
neutrinos varies significantly around the tracking regime, 
it has changed very little around today. This is due to 
the fact that the coupled system of neutrinos and dark 
energy has been stringently constrained by the current 
observations. It is trivial if the system has entered the 
tracking regime today when the fraction of the neutrino 
energy density is negligible |j2^. Meanwhile as the track- 
ing regime has been kept in the far future from today, the 
mass of neutrinos has changed little till the present epoch 
due to the slow rolling of (j). In this sense our scenario 
is hard to distinguish experimentally from the standard 
case with constant neutrino masses. On the other hand, 
however, if we do not expect a neutrino dominant future 
universe and the present epoch is closer to the track- 
ing regime, neutrino mass could have evolved from the 
past to now, as shown in the example of V = Vo<p^ and 
M = Yig. El We could then expect some 

distinguishable imprints from the uncoupled case|2lj. 

In addition, if the neutrino mass is produced through 
the seesaw mechanism and the scalar field only couples 
with the right-handed neutrinos, this scenario will change 
the predictions of the leptogenesis mechanism. As shown 

^ For the phantom scalar field, the effective Lagrangian after in- 
cluding the quantum effects from the derivative couplings is given 

by 

^./7 = ^ (-l + YiX^)aM09''</'. (26) 



explicitly in Ref. 10] the reheating temperature required 
by leptogenesis can be lowered and the bound on the light 
neutrino mass set by leptogenesis can also be relaxed. 

Finally, the derivative coupling in Eq. Q contributes 
also to the kinetic term. It is easy to show that the 
quantum effect at zero temperature vanishes. However, 
at finite temperature T the corrections are given by 

^correct = ^ (y^) ^^^^'^^ ■ (25) 

The effect is negligible at low temperature T <C A. How- 
ever, it could be important at very high temperatures 
when T ~ A. Ref.[lJ| has also considered another type 
of the back reaction caused by the derivative couplings 
in Eq. and shown explicitly the negligible effects at 
low temperature ^. 

We should point out that the main purpose of the cur- 
rent paper is the phenomenological study of the coupled 
system of dark energy and neutrinos. We do not try to 
build the dynamical relation of neutrinos and dark en- 
ergy and a possibility of neutrino condensate to achieve 
dark energy is pointed out in the footnote 1. Attempts 
for a possible dynamical link between the values of the 
neutrino mass and dark energy today have been studied 
in Refs. 5, 24]. Furthermore as shown in the previous 
studies of tracking dark energy models one usually needs 
some fine tuning for the energy density of quintessence 
today 0. 

In summary we have presented a systematic study on 
the cosmological evolution of the dark energy model with 
a coupled system including a dynamical scalar field (the 
quintessence) and the neutrinos. We have shown that the 
dynamics of this system drives the universe accelerating 
at present and leads to various interesting possibilities on 
the evolution of the universe. In our models, the neutri- 
nos may dominate the universe in the future. However, in 
most scenarios, the neutrinos will eventually decay away. 
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At high temperatures the system is a canonical scalar field while 
at low temperature a 'phase transition' occurs and the system 
turns to the phase of phantom. This tr ansi tion is somewhat 
favored by the present SN la observations \2?\ . 
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